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Abstract 

The one-dimensional XXZ model is studied in the presence of disorder in 
the Heisenberg Exchange Integral. Recent predictions obtained from renor- 
malization group calculations are investigated numerically using a Lanczos 
algorithm on chains of up to 18 sites. It is found that in the presence of 
strong X-Y-symmetric random exchange couplings, a "random singlet" phase 
with quasi-long-range order in the spin-spin correlations persists. As the pla- 
nar anisotropy is varied, the full zero-temperature phase diagram is obtained 
and compared with predictions of Doty and Fisher [Phys. Rev. B 45 , 2167 
(1992)]. 
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The study of quantum models in the presence of disorder is an emerging field onto which 
much attention has been focussed lately. Since all experimentally accessible systemsS are 
to some extent affected by randomness in the form of impurities, fields, or couplings, a 
thorough understanding of disorder effects can help in comparing experimental observations 
and theoretical predictions. In particular, weakly disordered, low-dimensional quantum 
spin systems are of interest, since the interplay between randomness and strong quantum 
fluctuations can be observed.! At T=0, phase transitions in quantum spin models are driven 
by zero-point fluctuations, as opposed to thermal excitations in their classical counterparts. 
However, when a random potential is introduced, phase transitions can be also driven by 
random fluctuations. This mechanism is particularly interesting in the case of marginally 
ordered systems, where the long-range Neel order in the 2D isotropic Heisenberg model has 
been found to be unstable towards thermal fluctuations and random fields, but not towards 
randomness in the exchange couplings.! 

The anisotropic spin-1/2 Antiferromagnetic Heisenberg Chain is a generic model of 
strongly correlated electrons. It is described by the Hamiltonian, 



where the notation is standard. Due to the low dimensionality, quantum fluctuations destroy 
long-range order in the region — 1 < A < 1, and the spin-spin correlations decay spatially 
with a power-law. Beyond the Heisenberg point (i.e. A > 1), a gap opens in the excitation 
spectrum and the system develops long-range Neel order with exponentially decaying corre- 
lation functions, while for A < — 1 there is a ferromagnetic region with Ising-type long-range 
order. 

Let us now introduce disorder in the form of X-Y symmetric random exchange couplings, 
i.e. such that the planar symmetry of H is not broken by the random potential, 



The random couplings Si are drawn from a uniform distribution P(5i) = 9 (Si — S J xy )8(S J ' xy — 
Si), where {Si) = and {(Si) 2 ) = 2(SJ xy ) 3 /3. The cut-off parameter SJ xy serves as a measure 
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for the strength of the random potential. The physical properties induced by this distribution 
are believed to be universal. However, in order to test this idea, we also studied random 



Here, a xy serves as a measure of the random strength. 

The properties of XXZ chains in the presence of various random potentials have recently 
been studied by C. A. Doty and D. S. Fisher using renormalization group techniques.^ It 
was found that, while random transverse fields destroy the (quasi) -long-range spin order, a 
power-law decay of the spin correlations may persist in the presence of random exchange 
couplings as long as the random Hamiltonian does not break the planar symmetry of H . 
In particular, it was predicted that a quasi-long-range-ordered phase extends from the X-Y 
regime (—1 < A < 1), when H random is switched on. 

In our study of the above system, we numerically diagonalized chains of up to 18 sites with 
periodic boundary conditions using a Lanczos algorithm. The observables were obtained 
from a quenched average, i.e. the ground state |0o(j)) of a chain was obtained for a given set 
of random couplings j = {Si}, and then the expectation value of some particular operators O 
were studied. This procedure was repeated for m ~ 500 different sets of random couplings, 
and finally the algebraic average over all m random samples was taken. The quenched 
average of an operator O is thus defined by 



"° 3=1 

First, we would like to address the question of whether quasi- long-range order persists 
in the region — 1 < A < 1 when the disorder potential H ran d om is switched on. The relevant 
observable is the normalized real-space spin-spin correlation function 

where N denotes the number of sites, and S = 1/2 in our study. 

In Fig.l, the spin-spin correlations u z (N/2) at the maximum separation (7 = N/2) are 
plotted as a function of the lattice size N at planar anisotropy A = 0.5 for a couple of 



exchange couplings drawn from a Gaussian distribution, P(5{) 



i 



exp (-5f/2a xy ). 




(3) 
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random strengths 5J xy . If the correlations decay with a power-law |co> z (/)| oc l~ Vz , we expect 
a straight line with negative slope r] z in a double-logarithmic plot. It is found that for all 
random strengths, 5J xy , a power-law decay (solid line) fits the numerical data much better 
than an exponential decay (dashed line), e.g. the x 2 ~ value obtained from least-square fits is 
typically two orders of magnitude larger when an exponential decay |co> z (/)| oc exp (— £1) is 
assumed. We observed a similar power-law behavior in a large region of parameter space. 

Why does the random potential not destroy quasi-long-range order in this region? Ac- 
cording to Doty and Fisher the "random singlet" phase which extends from the X-Y phase 
of the pure system (Hq) can be pictured in terms of randomly distributed tightly coupled 
singlet pairs of spins.! Those spins which are not bound in a singlet pair interact via virtual 
excitations. It turns out that these "almost-free" spins are anomalously strongly corre- 
lated. The probability that "almost-free" spins separated by a distance R interact strongly 
is proportional to 1/R 2 . This gives rise to the observed power-law behavior in the spin-spin 
correlations. The decay exponent is found to be r] z = 2. Note that also, in the exactly 
solvable X-Y limit (A = 0) the system maps into a tight-binding model of free fermions with 
random nearest-neighbor hopping. In this limit the decay exponent is given by rj z = 2 if a 
single characteristic localization length is assumed for the properties of the low-energy wave 
functions.^ 

In Fig. 2, we show rj z obtained in our numerical analysis, as a function of the disorder 
parameter SJ xy for various anisotropies A. The exponent has been extracted using chains 
of size N=6, 10, 14 and 18.1 The inset of Fig. 2 (a) shows the decay exponent r] z for the 
pure system Hq as it has also been obtained in Ref. 6. The exact diagonalization results 
are in excellent agreement with predictions from conformal invariance,! and in particular 
the Heisenberg limit (|u; z (/)| oc Z -1 ) and the X-Y limit (|u; z (/)| oc l~ 2 ) are nicely recovered. 
For negative anisotropies, (—1 < A < 0) conformal invariance predicts a constant exponent 
r] z = 2, which is also in reasonable agreement with our data, showing that our techniques 
can reproduce known results very accurately. 

On our finite chains and as we depart from the 5 J xy = limit, three regions can be 
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identified: 

(1) In the regime of small randomness (5J xy < J) the exponent rj z increases slightly as a 
function of the disorder parameter SJ xy , which is a sign of reduced order. 

(2) Around S J xy = J, there is an area of high competition between the quantum fluctua- 
tions of the original Hamiltonian (J J2i(S* Sf +1 + SfSf +1 ) and H random . Locally the random 
terms can compensate the zero-point fluctuations leading to an antiferromagnetic Ising-like 
behavior in the correlation functions. As a result, the decay exponent r] z has a dent with 
onset at around SJ xy = J, indicating a crossover into a more ordered Ising-like regime, where 
correlations decay more slowly than for the uniform system. 

(3) For large disorder, (8J xy » J) H random is the dominant term. The dependence of 
the decay exponent on the planar anisotropy in H becomes negligible, and it approaches 
f] z = 2 for all values of A, as it has been predicted by renormalization group arguments.^ 

In the vicinity of A = —0.75 the exponent r\ z behaves anomalously for small disorder. 
The observed decay in r] z for 8J xy between J and 2 J is due to ferromagnetic behavior in 
the real space spin-spin correlations. This anomaly is observed specially for anisotropies 
— 1 < A < —0.5. The dent of r\ z around 5J xy = J can be understood as a crossover into a 
phase of higher order. In particular, for A = —0.75 we observed a transition into a partially 
polarized phase indicated by the change of sign in the energy difference SE = E(S^ ot = 
0) — E(St ot = 1), where E(S^ ot = n) is the quenched ground state energy in the subspace 
with Sl ot = n. The inset of Fig. 2(b) shows SE as a function of the disorder parameter SJ xy 
at anisotropy A = —0.75 for a 14-site chain. It can be nicely seen that the transition into 
the partially polarized phase (0.55J < SJ xy < 3.05J) corresponds to the dent in r\ z in the 
same regime of disorder. 

In Fig. 3 (a), the dependence of the energy on the disorder parameters SJ xy and o xy 
at various anisotropies is shown for a 14-site chain. As the random potential becomes 
dominant, the system is allowed to relax into a ground state of higher entropy. The ground 
state energy drops proportionally to SJ xy (u xy ) in this region. In Fig. 3(b), we show how 
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the static structure factor (S zz (k) = J2j exp(— ikj)(0\S z S z +1 \0)) behaves as a function of the 
disorder parameters at antiferromagnetic momentum transfer k = 7r for the 14-site chain. 
In analogy to Fig.2, three regions can be identified. At low disorder the structure factor 
remains approximately unchanged. In the region of competition, Neel order is favored for 
positive anisotropies (0 < A < 1), resulting in an increase of the antiferromagnetic structure 
factor especially in the vicinity of A ~ 1. For negative anisotropies (— 1 < A < 0), the ditch 
in S zz (tt) indicates a crossover into a ferromagnetically polarized region. For large disorder, 
S zz (tt) becomes independent of A, and approaches the X-Y limit for all anisotropies. 

The boundary between the long-range-ordered regime and the "random singlet" phase 
is obtained from the correlations u z (N/2). In the "random singlet" phase, the spin-spin 
correlations at distance N vanish in the bulk limit as N — > oo. However, as the anisotropy is 
tuned across the critical value A c , u z (N/2) becomes finite, approaching |a; z (iVy2)| = 1 in the 
extreme Ising limit (A = oo). At zero disorder the Heisenberg point A c = 1 is nicely recovered 
as the critical point (Fig.3(c)). In Fig.3(d), we see that the transition point between these 
two phases is reduced to about A c = 0.75 at SJ xy = J.0 As a result of the strong competition 
effects in the region 5J xy ~ J, the antiferromagnetic phase bends into the random singlet 
regime in a "reentrant" transition, indicating a stronger antiferromagnetic order in this 
region. The whole boundary between "random singlet" and Neel phase is plotted in the 
phase diagram given in Fig.4. 

Both the "random singlet" and the Neel phase lie in the S z otal = subspace. On the 
other hand, as the ferromagnetic limit is approached, there is a transition into a partially 
polarized phase, i.e. the ground state no longer has S^ otal = 0. This phase boundary, as 
well as the transition from the partially into the fully polarized regime, is extracted from 
comparing the lowest energies of the various S^ otal subspaces (averaged over the ensemble 
of random couplings). In the region of competition between quantum fluctuations and the 
disorder term, the partially polarized phase bends into the "random singlet" regime, in 
analogy to the effect at the phase boundary between the "random singlet" and the Neel 
phase, as shown in Fig.4. 
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For low disorder, our results agree qualitatively with those of Doty and Fisher.Q However, 
their study predicts an X-Y-like "mole hill" phase in the region — 1 < A < —0.5, and for small 
disorder. Numerically, it is hard to distinguish this "mole-hill" from the "random singlet" 
regime, because both phases show power-law behavior in the correlation functions. However, 
from our exact diagonalization data we have observed a region (denoted with a question mark 
in Fig. 4) which has power-law decay, and is a member of the S* otal = subspace, but does 
not have any remnant antiferromagnetic correlations, as has been discussed above in the 
inset of Fig.2(b). We are currently investigating, whether this regime can be identified with 
the "mole hill" predicted by Doty and Fisher. 

In summary, we have presented the first numerical study of the spin-1/2 XXZ chain in 
the presence of a random exchange potential {H random ). In contrast to a random fieldjl 
quasi- long-range order of the zero-disorder X-Y regime — 1 < A < 1 is not destroyed by an 
X-Y symmetric random exchange. Also, Ising-type long-range order persists in the presence 
of small random exchange couplings. The power-law behavior in the "random singlet" phase 
may be due to virtual interactions of "almost-free" spins which are not bound in randomly 
strong singlet pairs. A complete phase diagram is provided. In addition, we have found 
an interesting reentrant transition of the ordered phases (in both the ferromagnetic and 
antiferromagnetic Heisenberg limits) when exchange disorder is included. Such a novel type 
of behavior (order induced by random couplings) deserves additional study. 
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FIGURES 

FIG. 1. : Double-logarithmic plot of real-space spin-spin correlations |^ Z (Z)| at maximum sep- 
aration (I = N/2) as a function of lattice size. The squares represent data obtained from exact 
diagonalizations, the solid lines are fits to power-law decay |^ Z (0I = Al~ Vz and the dashed lines 
are fits to an exponential decay |w z (Z)| = Aexp (—£1). The size of the squares is comparable to the 
magnitude of the corresponding error bars. 

FIG. 2. : (a) Exponents of the power-law decay |<-J 2 (0I = Al~ riz as a function of the disorder 
parameter SJ xy for various positive planar anisotropics. The inset shows 7] z as a function of 
anisotropy in the limit of no disorder, (b) Same as (a) but for negative anisotropics. The inset 
shows the energy difference 5E = E(S^ ot = 0) — E(Sf ot = 1) as a function of SJ xy for the 14-site 
chain. The change in the sign of 5E indicates the presence of a partially polarized phase for 
0.55J < SJ xy < 3.05 J at anisotropy A = —0.75. 

FIG. 3. : (a) Ground state energy of the 14-site spin-1/2 XXZ chain as a function of the 
disorder parameter 8J xy at various planar anisotropies. The random exchange couplings are drawn 
from a uniform distribution with cut-off 5J xy . The inset shows the same except when the random 
exchange couplings are obtained from a Gaussian distribution of width a xy . (b) Antiferromagnetic 
structure factor vs. disorder parameter 5J xy for the 14-site spin-1/2 XXZ chain. The inset shows 
the same but when the random exchange couplings are obtained from a Gaussian distribution of 
width a xy . (c) Real space correlation functions at the maximum distance for an N site chain as 
a function of anisotropy at zero disorder. The bulk limit N = oo is extracted from a finite size 
study, (d) Same as (c) at disorder 5J xy = J. 

FIG. 4. : Phase diagram of the spin-1/2 XXZ chain in the presence of a random exchange 
potential. The question mark denotes the "mole hill" phase discussed in the text. 
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